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Abstract. Measuring the deviation from scale invariance of the primordial power spectrum 
is a critical test of inflation. In this paper we reconstruct the shape of the primordial power 
spectrum of curvature perturbations from the cosmic microwave background data, including 
the 7-year Wilkinson Microwave Anisotropy Probe data and the Atacama Cosmology Tele- 
scope 148 GHz data, by using a binning method of a cubic spline interpolation in log-log 
space. We find that the power-law spectrum is preferred by the data and that the Harrison- 
Zel'dovich spectrum is disfavored at 95% confidence level. These conclusions hold with and 
without allowing for tensor modes, however the simpler model without tensors is preferred 
by the data. We do not find evidence for a feature in the primordial power spectrum — in 
full agreement with generic predictions from cosmological inflation. 
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1 Introduction 


Inflation in the early Universe predicts primordial power spectra of curvature and tensor 
perturbations which lead to a distinct imprint on the observed angular power spectrum of 
the cosmic microwave background (CMB) radiation. In the simplest models of slow-roll 
inflation the power spectra can analytically be expressed in terms of the inflaton potential 
and its first derivative at the Hubble scale crossing during inflation. The angular power 
spectrum of the CMB is a convolution of the primordial power spectrum with a radiative 
transfer function, which can be obtained by integrating the full Einstein-Boltzmann system of 
differential equations. Therefore, given an inflaton potential or a primordial power spectrum, 
the angular power spectrum of the CMB can numerically be calculated, and thus the inflaton 
potential or the primordial power spectrum can be fitted to CMB data [1-3]. 

The simplest parameterization of the primordial power spectrum is the phenomeno- 
logical Harrison-Zel’dovich (scale invariant) spectrum with one free parameter. The next 
simplest ansatz is to assume that there is no distinguished scale, which leads to a simple 
power-law parameterization. For this power-law ansatz, the scale-invariant spectrum is ex- 
cluded at more than 99% confidence level by the Wilkinson Microwave Anisotropy Probe 
(WMAP) data [2]. Moreover, a slightly tilted power-law spectrum is still an excellent fit to 
the data even adding a running index which characterizes the deviation from a pure power 
law [2]. Other specific parameterizations of the primordial power spectrum, motivated by 
theoretical models or features of the observed data, have been considered: for example, a 
broken power spectrum [4] due perhaps to an interruption of the inflaton potential [5], and 
a cutoff at large scales [6, 7] motivated by suppression of the lower multipoles in the CMB 
anisotropies [8, 9]. 

Since a modest deviation from scale invariance of the primordial power spectrum is a 
critical prediction of inflationary models, it is a vital test of the inflationary paradigm. Exact 
scale invariance would be expected for a perfectly stationary epoch of inflation, which would 
never end. On the other hand, a strong deviation from scale invariance could falsify the 
idea of cosmological inflation. Thus it is important to probe the shape of the primordial 
power spectrum. A strong theory prior on it could lead to misinterpretation and biases in 
parameter determination. Given our complete ignorance of the underlying physics of the 
very early Universe, some model-independent approaches to reconstruct the shape of the 
primordial power spectrum from existing data have been employed, involving for example 
linear interpolation [10], a minimally-parametric reconstruction [11], wavelet expansions [12], 
principle component analysis [13], and a direct reconstruction via deconvolution methods [14— 
16]. Clearly, multiple methods are needed to cross-check each other and to contribute their 
respective strengths to our understanding. 
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In this work, we reconstruct the shape of the primordial power spectrum using a bin- 
ning method of a cubic spline interpolation in the logarithmic amplitude and logarithmic 
wavenumber space, where the power-law spectrum as a special case is just a straight line. 
This method guarantees the positivity of the power spectrum. We use not only the seven- 
year WMAP data [2] but also small-scale CMB data from the Atacama Cosmology Telescope 
(ACT) experiment [3], with two main astrophysical priors on the Hubble constant (Ho) mea- 
sured from the magnitude-redshift relation of 240 low-z Type Ia supernovae at z < 0.1 [17] 
and on the distance ratios of the comoving sound horizon to the angular diameter distances 
from the Baryon Acoustic Oscillation (BAO) in the distribution of galaxies [18]. Using data 
from ACT besides WMAP allows us to extend the measured k range to smaller scales. The 
combination of both data sets provides significant improvements in the reconstruction of the 
spectrum. 

In previous works, tensor modes have been commonly ignored when testing for the 
shape of the primordial power spectrum. Here we consider the contribution from tensor 
modes and assume a scale-invariant spectrum of tensor perturbations. We compare different 
spectral shapes by means of the Akaike information criterion, including a comparison between 
models with and without tensor perturbations. We show that current data prefer a power-law 
shape and no tensor perturbations. 

This paper is organized as follows. In Section 2, we describe the method and the data 
used in this analysis. In Section 3, we present our results. Section 4 is devoted to conclusions. 


2 Analysis method 


We divide the primordial power spectrum of curvature perturbations, Pr, into Npin bins 
equally spaced in Ink between kı = 0.0002 Mpc! and knpa = 0.2 Mpc™t for the 7-year 
WMAP data with l < 1200 for the TT power spectrum and / < 800 for the TE power spec- 
trum. Moreover, we also use the low-/ temperature and polarization angular power spectra 
from WMAP. Adding the ACT data between 1000 < | < 3000 allows a wide wavenumber 
range with ky,,, = 0.3 Mpc™t. 

To reconstruct a smooth spectrum with continuous first and second derivatives with 
respect to Ink, we use a cubic spline interpolation to determine logarithmic values of the 
primordial power spectrum for Ink; < Ink < Inkj,,. Here boundary conditions are adopted, 
where the second derivative is set to zero. For k < kı or k > kyy,,,, we fix the slope of the 
primordial power spectrum at the boundaries since the CMB data place only weak constraints 
on them. This reconstruction method has three advantages: firstly, it is easy to detect 
deviations from a scale-invariant or a power-law spectrum because both the scale-invariant 
and power-law spectra are just straight lines in the In k-In Pr plane. Secondly, negative values 
of the spectrum can be avoided by using ln Pr(k) instead of Pr(k) for splines with steep 
slopes. Thirdly, the shape of the power spectrum reduces to the scale-invariant or power-law 
spectrum as a special case when Np, = 1,2, respectively. To summarize, we consider the 
following form of the power spectrum in this work, 


a | In È + In Pg(ki), k< ki; 
1 
In Pr(ki), k € {ki}; 
In Pr(k) = cubic spline, kel kx kisis (2.1) 
dln Pr (k 
Aa ) i ln TA +n PR EN yia) KE > ENyn: 


bin 
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We consider a spatially flat ACDM Universe described by Npin primordial spectrum 
parameters A; = ln [100 Pp (ki)] and four background parameters (Q)h?, Q-h?, Os, T), 
where Q,h? and Q,h? are the physical baryon and cold dark matter densities relative to the 
critical density, ©, is the ratio of the sound horizon to the angular diameter distance at 
decoupling, and 7 is the reionization optical depth. We also consider the Sunyaev-Zel’dovich 
(SZ) effect, in which CMB photons scatter off hot electrons in clusters of galaxies. Given a 
SZ template it is described by a SZ template amplitude Asz as in the WMAP papers [1, 2]. 
For the 148 GHz ACT data [3], aside from Agz we use, following the ACT analysis, two more 
secondary parameters: A, and Ae. The former is the total Poisson power at | = 3000 from 
radio and infrared point sources. The latter is the template amplitude of the clustered power 
from infrared point sources. If a contribution to the CMB from tensor modes is considered, 
we assume a scale-invariant primordial power spectrum of tensor perturbations. Current 
CMB measurements are insensitive to deviations from this choice, as the tensor modes are 
subdominant in the TT angular power spectrum and contribute at most up to l ~ 100. 

Our analysis is carried out using a modified version of the publicly available CosmoMC 
package, which explores the parameter space by means of Monte Carlo Markov Chains [19]. 
Besides the 7-year WMAP data (J < 1200) including the low-l temperature (2 < l < 32) 
and polarization (2 < l < 23) data, we use two main astrophysical priors: the present-day 
Hubble constant Ho from the magnitude-redshift relation of 240 low-z Type Ia supernovae 
at z < 0.1 [17], and the angular diameter distances out to z = 0.2 and 0.35, measured from 
the two-degree field galaxy redshift survey and the sloan digital sky survey data [18]. To 
enlarge the range in k, we use the WMAP data in combination with the 148 GHz ACT data 
during its 2008 season. For the ACT data, we focus on the band powers in the multiple 
range 1000 < l < 3000. Following Ref. [3] for computational efficiency the CMB is set to 
zero above | = 4000 where the contribution is subdominant, less than 5% of the total power. 
To use the 148 GHz ACT likelihood there are three secondary parameters, Asz, Ap and Ae. 
We impose positivity priors on these parameters and use the SZ template and the clustered 
source template provided by the ACT likelihood package. 

To judge whether adding more bins effectively improves the fit, we apply the Akaike 
information criterion (AIC) defined by 


AIG = —21n Lmax + 2Npar, (2.2) 


where Lmax is the maximum likelihood achievable by the model and Npar the number of 
parameters. This criterion arises from an approximation to the Kullback-Leibler information 
entropy, used to measure the difference in goodness of fit between two models. It sets up a 
tension between goodness of fit and complexity of the model, with the best model minimizing 
the AIC. Typically, models with too few parameters give a poor fit to the data and hence have 
a low log-likelihood, while those with too many are penalized by the second term. Although 
it is not always clear how big a difference in AIC is required for the worse model to be 
significantly disfavored, typically a difference of 6 or more should be taken seriously. This 
information criterion was used to carry out cosmological model selection [20]. 


3 Results 


Fig. 1 shows the reconstructed power spectrum of curvature perturbations without tensor 
modes, derived from WMAP7+Ho+BAO (top panels) and WMAP7+ACT+Ho+BAO (mid- 
dle and bottom panels). Large uncertainties of the reconstructed spectrum at low-l mainly 
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arise from the cosmic variance while at high-/ different systematics affect the reconstructed 
spectrum due to point source emission and the SZ effect. As the number of bins increases, 
the error bars of individual bins grow due to the increase in degrees of freedom. We can see 
that the power-law spectrum is an excellent fit to the data while the scale-invariant spectrum 
is excluded by 95% confidence level at small scales if tensor modes are ignored. Including the 
ACT data allows us to increase kN,» to 0.3 Mpc7!. At this scale the scale-invariant spectrum 
is far outside of the 2o error bound. As shown in Table 1, with the increasing number of 
bins the AIC value first decreases and then increases for Npin > 2. We find that Npin = 2 
(a power-law shape) is the best model with and without ACT data and Npin = 1 (Harrison- 
Zel'dovich) provides the least favored fit to the data. To demonstrate strong constraints on 
the spectrum at small scales from the ACT data, we also reconstruct the spectrum in the 
case of kN, = 0.2 Mpc~? in the middle panels of Fig. 1. Compared to the top panels, the 
errors in the last bin are reduced by including the ACT data. 
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Figure 1. Reconstruction of the primordial power spectrum without tensor modes using a cubic spline 
interpolation for different binnings, derived from the WMAP7+Ho+BAO combination (top panels) 
and from the WMAP7+ACT+H)+BAO combination for two different choices of ky,,,, (middle and 
bottom panels). The dotted and dashed lines represent the scale-invariant (Npin = 1) and power- 
law (Nbin = 2) spectra respectively. The blue and red error bars show lo and 20 uncertainties 
respectively. Note also that the best fit optical depth depends quite strongly on the assumed shape 
of the primordial spectrum. 
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Data Nbin Npar — ln Lmax Ax? AIC AAIC 
WMAP7+ Hp +BAO 1 6 3742.4 0 7496.8 0 
2 7 3738.7 —74 7491.4 —5.4 
3 8 3738.4 —8.0 7492.8 —4.0 
4 9 3738.4 —8.0 7494.8 —2.0 
5 10 3738.4 —8.0 7496.8 0 
+ACT 1 8 3752.1 0 7520.2 0 
2 9 3747.9 —8.4 7513.8 —6.4 
3 10 3747.1 —10.0 7514.2 —6.0 
4 11 3747.2 —9.8 75164 —3.8 
5 12 3747.2 —9.8 75184 —1.8 


Table 1. The y?, AIC values and their differences with respect to the scale-invariant power spectrum 
for various models, derived from WMAP7+H)+BAO and from WMAP7+ACT+H)+BAO. Here the 
contribution from tensor modes are ignored. 


Assuming a scale-invariant spectrum of tensor perturbations, the reconstructed primor- 
dial power spectrum of curvature perturbations are shown in Fig. 2. Due to the contribution 
of tensor modes to the angular power spectra of the CMB, the amplitude of the Harrison- 
Zel'dovich spectrum is smaller than one in Fig. 1. For WMAP7+AHo+BAO, there is no 
convincing deviation from either a simple scale-invariant or a power-law spectrum as shown 
in the top panels of Fig. 2. As expected, including the ACT data will improve the measure- 
ments of the spectrum. In the case of Npin = 3 the reconstructed spectrum deviates from 
the scale-invariant one on small scales at 95% confidence level. In the Npin = 4 case, since 
the error bars become larger, it becomes consistent with the scale-invariant spectrum. From 
Table 2, we can see that increasing the number of bins beyond 2 increases the AIC value 
slightly. Hence, under the assumption of a scale-invariant tensor spectrum, the power-law 
spectrum provides a better fit to the data. Comparing Table 2 with Table 1, we find that the 
scale-invariant spectrum with tensor modes is more favored than one without tensor modes 
while the power-law spectrum without tensor modes is more favored than one with tensor 
modes. 

In Table 3 we summarize the estimated primordial power spectrum values as well as 
the other cosmological parameters. Note that the best fit reionization optical depth depends 
quite strongly on the assumed shape of the primordial spectrum. As Npin increases from 1 
to 3, the reionization optical depth first decreases and then increases drastically. 


4 Conclusions 


In this work we have reconstructed the smooth shape of the primordial power spectrum of 
curvature perturbations without or with tensor modes from the 7-year WMAP data in com- 
bination with the small-scale CMB data from the ACT experiment. We adopt the cubic 
spline interpolation method in the In k-ln Pr plane to avoid negative values of the ampli- 
tude and to make the detection of deviations from the power law easy. If we ignore tensor 
modes, the scale-invariant spectrum is excluded by the combined data at 95% confidence 
level, which is consistent with the result of Ref. [2]. Moreover, we find no convincing de- 
viation from a simple power-law spectrum as suggested in Refs. [11, 15]. If tensor modes 
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Figure 2. Reconstruction of the primordial curvature power spectrum with a scale-invariant spec- 
trum of tensor perturbations using a cubic spline interpolation for different binnings, derived from 
WMAP7+Ho+BAO (top panels) and from WMAP7+ACT+Ho+BAO for two different choices of 
kn,;, (middle and bottom panels). The dotted and dashed lines represent the scale-invariant and 
power-law spectra respectively. The blue and red error bars show lo and 20 uncertainties respec- 
tively. 


are included, for WMAP7+H)+BAO the scale-invariant spectrum lies within the 20 bound. 
For WMAP7+ACT-+ Ho+BAO, there a deviation from the scale-invariant spectrum at small 
scales is found at the 95% confidence level under the assumption of a scale-invariant tensor 
spectrum. 

Comparing Fig. 1 with Fig. 2, we find that including tensor modes enhances the bending 
of the primordial power spectrum. Moreover, the AIC values listed in Table 1 and Table 2 
indicate that the power-law spectrum without tensor modes is the best fit to the combined 
data. We do not find evidence for a feature in the primordial power spectra, although a bent 
spectrum with a suppression of low k and high k modes with respect to a power-law spectrum 
is only slightly disfavored by the AIC (—6.4 vs. —6.0 without tensors and —4.4 vs. —4.0 with 
tensors). Without ACT the case for a bend is weaker. 

We emphasize that our reconstruction method is insensitive to local features in the 
primordial power spectrum, but is sensitive to the overall shape, since the cubic spline in- 
terpolation makes the primordial power spectrum smoother than linear interpolation [10]. 
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Data Nbin Npar — ln Lmax Ax? AIC AAIC 

WMAP7+ Hp +BAO 1 7 3740.5 —3.8 7495.0 —1.8 
2 8 3738.8 —7.2 7493.6 —3.2 
3 9 3738.4 —8.0 7494.8 —2.0 
4 10 3738.4 —8.0 7496.8 0 
5 11 3738.2 —8.4 7498.4 1.6 

+ACT 1 9 3750.4 —3.4 7518.8 —1.4 
2 10 3747.9 —8.4 7515.8 —4.4 
3 11 3747.1 —10.0 7516.2 —4.0 
4 12 3747.1 —10.0 7518.2 —2.0 
5 13 3747.2 —9.8 7520.4 0.2 


Table 2. The y’, AIC values and their differences with respect to the scale-invariant power spectrum 
without tensor modes for various models with tensor modes, derived from WMAP7+H)+BAO and 
from WMAP7+ACT+H,+BAO. We assume a scale-invariant spectrum of tensor perturbations. 


Nbin Aj A2 A3 A4 As Qah Qeh Qa T 
WMAP7+H0o+BAO (ENpin = 0.2 Mpc!) 

1 3.132+0.031 = = = = 0.02374 0.1131 0.738 0.101 
2 3.269+0.060 3.037+0.045 — — — 0.02261 0.1125 0.727 0.088 
3 3.096+0.167 3.168+0.032 3.005+0.054 — — 0.02224 0.1137 0.719 0.092 
4 2.990+0.296 3.172+0.042 3.132+0.038 3.031+0.076 = 0.02233 0.1130 0.722 0.094 
5 3.057+0.505 3.137+0.076 3.176+0.039 3.105+0.042 3.054+0.099 0.02224 0.1127 0.722 0.095 
+Tensor 

1 3.109+0.033 = = = = 0.02364 0.1117 0.741 0.094 
2 3.211+0.074  3.052+0.048 = = = 0.02289 0.1117 0.733 0.087 
3 2.822+0.267 3.145+0.038  3.009+0.055 — — 0.02244 0.1137 0.722 0.097 
4 2.772+0.380 3.084+0.077 3.141+0.041 3.017+0.079 — 0.02251 0.1132 0.724 0.098 
5 3.226+0.549 2.837+0.227 3.115+0.057 3.102+0.042 3.076+0.100 0.02219 0.1124 0.724 0.095 
WMAP7+ACT+Ho0+BAO (E Nyin =0.3 Mpc‘) 

1 3.133+0.032 — — — — 0.02358 0.1124 0.741 0.102 
2 3.279+0.056  3.014+0.047 = = = 0.02239 0.1125 0.727 0.085 
3 3.072+0.148 3.173+0.036 2.967+0.056 = = 0.02203 0.1143 0.716 0.094 
4 3.033+0.274  3.168+0.039 3.136+0.037  2.965+0.080 — 0.02204 0.1146 0.715 0.094 
5 3.028+0.452 3.164+0.068 3.176+0.039 3.104+0.038 2.957+0.099 0.02210 0.1148 0.714 0.095 
+Tensor 

iL 3.110+0.031 — — — — 0.02347 0.1109 0.745 0.095 
2 3.241+0.048 3.034+0.037 — — — 0.02253 0.1117 0.731 0.086 
3 2.800+0.258 3.158+0.039 2.961+0.062 = = 0.02225 0.1146 0.718 0.098 
4 2.793+0.374 3.098+0.069 3.146+0.039 2.956+0.080 = 0.02227 0.1148 0.717 0.098 
5 3.080+0.498 2.941+0.173 3.142+0.044 3.102+0.039 2.980+0.102 0.02204 0.1152 0.713 0.094 


Table 3. Mean values and marginalized 68% confidence levels for the primordial spectrum parameters, 
and the mean values for the background parameters. 


Therefore, it is complementary to wavelet expansions [12] and principle component analy- 
sis [13], which provide rigorous methods to search for sharp features. This type of analysis 
is also complementary to the direct testing of slow-roll inflation [21, 22] and the analysis of 
the WMAP team [1, 2]. 
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